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ABSTRACT 

In  June  1988,  the  Naval  Ocean  Research  and  Development  Activity  (NORDA)  collected 
some  "in— water"  data  using  its  Towed  Underwater  Pumping  System  (TUPS)  in  the  near-shore 
waters  off  St.  Andrews  State  Park,  Shell  Island,  Florida.  These  in  situ  data  include  latitude; 
longitude;  depth  in  meters;  narrow-band  upwelling  at  465  nm,  507  nm,  and  532  nm;  broad-band 
downwelling  collected  at  the  surface;  temperature;  salinity;  and  transmissivity.  In  this  paper,  we 
investigate  the  relationship  between  depth  and  the  normalized  upwelling  irradiance  (upwelling 
divided  by  downwelling)  in  the  three  bands.  Algorithms  used  to  calculate  water  depth  from 
remotely  sensed  airborne  and  sate'lite  multispectral  data  are  applied  to  the  TUPS  data  and  results 
compared.  The  TUPS  data  have  the  advantage  over  most  aircraft—  and  satellite— collected  data 
because  they  were  collected  over  an  essentially  uniform  bottom  type  (smooth  sandy  bottom  with 
steady  slope)  and  have  no  atmospheric  contamination.  A  new  algorithm  for  depth  calculation  is 
proposed. 


1.  INTRODUCTION 

Over  the  past  few  years,  the  Naval  Oceanographic  and  Atmospheric  Research  Laboratory 
(NOARL),  formerly  the  Naval  Ocean  Research  and  Development  Activity  (NORDA),  has  been 
developing  and  testing  an  Airborne  Bathymetric  Survey  System  (ABS)  designed  to  measure  water 
depth  in  coastal  waters  (see  Hickman  et  al.5).  The  ABS  integrates  active  laser  soundings  with 
simultaneously  obtained,  passive,  multispectral  scanner  digital  imagery.  Algorithms  extrapolate 
laser  soundings  and  scanner  data  to  scene-wide  bathymetry.  If  sufficient  accuracies  can  be 
obtained,  then  this  technique  would  provide  for  inexpensive  and  rapid  production  of  bathymetric 
charts  for  shallow— water  areas. 

To  support  the  ABS  project,  to  provide  environmental  data  for  ground  truth,  and  to  help 
improve  ABS  algorithms,  NORDA  conducted  a  field  experiment  in  the  near-shore  waters  off  St. 
Andrews  State  Park  near  Panama  City,  Florida,  in  June  1988  (see  Guinasso  and  Wiesenburg4). 

In  this  field  experiment,  NORDA’s  Towed  Underwater  Pumping  System  (TUPS)  was  used  to 
collect  a  variety  of  environmental  data,  including  water  depth  and  upwelling  irradiance,  in  three 
narrow  spectral  bands  centered  at  465  nm,  507  nm,  and  532  nm.  Data  were  collected  on  the  south 
side  of  Shell  Island  (St.  Andrews  State  Park)  on  19  and  20  June  1988. 

In  this  paper  we  investigate  current  multispectral  algorithms  employed  for  shallow-water 
bathymetry  using  the  TUPS  data;  in  particular,  we  investigate  the  relationship  between  depth 
and  the  normalized  upwelling  irradiance  (upwelling  divided  by  downwelling)  in  the  three  bands. 
The  bathymetry  models  d;ccussed  here  are  based  on  the  fact  that  sunlight  is  exponentially 
attenuated  twice  as  it  passes  through  the  water  column  and  is  reflected  off  the  bottom.  Given  a 
data  base  of  known  depths  at  known  positions,  and  given  multispectral  imagery  that  encompasses 
the  data  base,  one  regresses  the  model  against  the  image  values  and  the  depth  data  base  to 
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determine  model  parameters,  which  then  can  be  applied  to  all  positions  in  the  imagery.  The 
TUPS  data  have  the  advantage  over  most  data  collected  by  aircraft  or  satellite  in  that  they  were 
collected  over  relatively  clear  water  with  an  essentially  uniform  bottom  type  (smooth  sandy 
bottom  with  steady  slope  from  the  beach  to  deeper  water  »  25  m)  and  have  no  atmospheric 
contamination.  Some  new  algorithms  for  depth  calculation  are  proposed  as  well.  We  also 
introduce  a  new  technique,  which  incorporates  an  effective  attenuation  coefficient  term  which 
varies  as  a  function  of  depth. 


2.  THE  TOWED  UNDERWATER  PUMPING  SYSTEM 

The  TUPS  contains  a  depth  sensor  (0  —  100  m  with  2— cm  resolution)  and  narrow  band 
(±10  nm)  upwelling  irradiance  sensors  at  wavelengths  465  nm,  507  nm,  and  532  mu.  NwRDA  s 
TUPS  Environmental  Sensor  System  (TESS)  was  used  to  collect  incident  light  measurements  for 
normalization  with  the  TUPS  data.  Employing  two  pyroheliometers,  broadband  downwelling 
irradiance  data  were  collected  simultaneously  with  the  TUPS  upwelling  data.  This  downwelling 
measurement  is  important,  since  clouds  overhead  reduce  the  amount  of  light  hitting  the  ocean 
surface  and  correspondingly  reduce  the  amount  of  upwelling  irradiance.  Table  1  shows  the  sensor 
configuration  of  TUPS,  and  Table  2  shows  the  TESS  configuration;  these  tables  are  adapted  from 
Guinasso  and  Wiesenburg4. 

Throughout  this  paper,  we  will  study  various  models  relating  normalized  upwelling  to 
water  depth.  The  light  signals  shall  be  normalized  irradiances  in  the  sense  that  we  take  = 

/pi,  where  ^  is  the  upwelling  radiance  (in  pW/cm2/nm)  received  on  the  TUPS  and  pi  is 

the  broadband  downwelling  radiance  measured  (in  /iW/cm2/nm)  at  the  surface  by  pyroheliometer 
1  in  the  TESS.  This  normalized  radiance  is  used  to  reduce  such  effects  as  clouds  passing  overhead. 
The  notation  LI,  L2,  L3  is  used  to  denote  these  normalized  upwelling  irradiances  in  wavelengths 
465  nm,  507  nm,  532  nm,  respectively.  Figure  1  shows  the  normalized  upwelling  versus  depth. 
Table  3  shows  part  of  the  correlation  matrix  of  TUPS  and  TESS  measurements;  for  the  full 
matrix  see  Guinasso  and  WiesenburgV  The  reader  should  note  that  the  two  pyroheliometers  are 
over  99%  correlated;  hence,  the  choice  of  which  to  use  for  the  normalization  is  immaterial.  Also 
note  that  the  three  nonnormalized  upwelling  bands  are  over  97%  correllated  with  each  other. 

Thus,  not  much  improvement  of  multiband  algorithms  over  single  band  algorithms  can  be 
expected. 


3.  THE  OPTICAL  BATHYMETRY  MODEL 

In  this  discussion,  we  assume  the  simple  radiance  transfer  model  based  on  the 
Lambert-Beers  Law,  called  herein  the  Single  Band  Model  (also  called  Bouguer’s  Law9  and  Jerlov’s 
Model7),  in  which  the  radiance  in  wavelength  A  at  water  depth  z  is  given  by  the  equation 

LA  =  LA„  +  RAe*P(-f-*A-*)-  (1) 

where 

L^  is  the  received  upwelling  radiance  value  in  the  wavelength  A, 

L^  is  the  observed  radiance  over  optically  deep  water, 

R^  is  a  radiance  term  that  is  sensitive  to  bottom  reflectance, 
is  an  effective  water  attenuation  coefficient  at  wavelength  A, 
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f  is  a  geometrical  constant  relating  the  angle  of  solar  incidence  and  the  sensor 
viewing  angle  and  f  =  2  for  nadir  view  with  sun  at  zenith. 

Solving  for  z  in  equation  (1)  yields 

z  =  (log(RA)  -  log(LA-LA J)/(f. *A)  (2) 

where  log(— )  stands  for  the  natural  logarithm.  The  difficulty  in  calculating  water  depth  directly 
from  equation  (2)  is  that  the  values  of  and  are  generally  not  known.  With  the  notation 

=  log(L^— L^),  the  depth  equation  is  recast  in  the  regression  form 

z  =  A  +  B-Xa,  (3) 

where  A  and  B  are  coefficients  to  be  determined  by  regressing  the  equation  against  known 
ground  truth,  minimizing  the  difference  between  the  known  depths  and  the  depths  calculated  from 
equation  (3)  using  measured  light  signals  and  L^.  In  order  to  do  this  regression,  it  is 

assumed  that  both  and  *A  are  constant  over  the  ground  truth. 

Paredes  and  Spero8  generalized  this  technique  into  a  multiband  algorithm  to  improve 
performance  by  reducing  the  sensitivity  of  the  algorithm  due  to  variation  in  bottom  reflection. 
They  suggested  that  it  is  possible  to  include  information  contributed  from  multiple  spectral  bands 
by  forming  a  weighted  sum  of  single-band  depth  solutions  as  given  by  equation  (3).  We  write  for 
n  wavelengths, 

n  n 

z  =  X  a,i'[(1/f,«i)'(l06(Ri)-Xi)].  with  ^0^=1.  (4) 

i  i 

This  immediately  generalizes  to  a  (regression)  model  in  the  form 


z  =  A0+XAi’Xi’  (5) 

i 

where  the  constants  A.  are  determined  to  minimize  the  error  between  the  model— predicted 

depths  and  the  measured  depths  at  the  ground  truth  control  points.  Clark  et  al.2  demonstrated 
this  model’s  improved  performance  using  Landsat  Thematic  Mapper  imagery,  and  experiments 
reported  in  Clark  et  al.3  add  further  evidence  to  this  improvement. 

Using  the  TUPS  normalized  upwelling  bands  LI,  L2,  and  L3,  we  performed  single-band 
regression  fits  of  equation  (3)  and  two-  and  three-band  regression  fits  of  equation  (4)  with  n  =  2 
and  3.  The  results  are  tabulated  in  Tables  4,  5,  and  6.  In  these  tables,  the  term  RM  stands  for 
the  statistical  (multiple)  correlation  coefficient  and  the  term  RMS  is  the  root  mean  square  of  the 
error  (calculated  depth  minus  calibration  depth). 

The  following  Figure  2  shows  scatter  plots  of  both  log(Ll)  and  log(Ll  -  Ll^).  Note 

that  the  scatter  plot  of  log(Ll)  has  roughly  constant  slope  for  3  £  z  £  10  m,  and  that  in  the 
scatter  plot  of  log(Ll  -  Llmin),  the  spread  of  the  plot  increases  dramatically  after  10  m.  This 
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suggests  that  «1  is  not  constant  and  that  the  errors  in  computing  the  bathymetry  will  increase  as 
we  take  the  computations  deeper.  Similar  effects  are  seen  in  scatter  plots  for  bands  L2  and  L3. 


4.  AN  EMPIRICAL  MODEL 

From  the  single— band  model,  a  graph  of  the  radiance  in  the  wavelength  A  versus  water 
depth  z  can  be  expected  to  have  the  appearance  illustrated  in  Figure  3.  Here,  we  assume  z 

is  sufficiently  deep  so  that  the  sign;  ’  received  at  this  depth  is  solely  due  to  scattering  within  the 
water  column  (no  bottom  reflectance).  Thus,  based  on  the  Lambert— Beers  Law,  we  expect  to  see 
a  curve  of  the  form 


where 


L^  =  A*exp(— 2*k^-z)  -F  B, 


Lim  L,  =  Lim  L,  =  L  •  =  B  and 
A  A  min 

z  “*  z  ^  z  ~>  oo 
max 

Lim  L,  =  L  --  A-exp(— 2- k, -z  .  )  +  L  ■  , 
A  max  FV  A  min'  min’ 

z  -*  z 


from  which  it  follows  that 

LA  =  Lmin  +  (Lmax-Lmin)-exp(-2-'tA-(z-zmin))’  <9) 

which  agrees  with  the  model  given  in  equation  (1)  by  setting  Lmin  =  L^  and  =  (Lmax  - 

Lmin)'exp(2'*A’zmin)-  Letting  AL  =  Lmax  “  Lmin  and  dropping  the  subscript  A,  we  have 
the  equation 

L  =  Lm,n  +  AL-e*P(-2-K-(z-zmin)>-  0°) 

which  we  shall  call  the  empirical  equation.  From  this  equation  it  is  easy  to  solve  for  k  obtaining, 
«=  0.5*[log(AL)-log(L-Lmin)]/(z-zmin).  (11) 

From  an  examination  of  the  values  of  LI,  we  empirically  obtained 

L1max  =  4  827  when  zmin  =  15  ra'  L1min  =  0  320  (12) 


Substituting  the  values  for  LI  into  equation  (11)  and  plotting  n  versus  depth  by  plugging  in  the 
TUPS  data  values,  we  obtained  Figure  4.  It  is  important  to  observe  that  here  nl  represents  an 
overall  effective  attenuation  coefficient  and  not  the  diffuse  attenuation  coefficient.  Moreover, 
these  empirical  «1  values  are  calculated  from  different  spatial  locations  that  have  various 
associated  depths  as  opposed  to  most  diffuse  attenuation  calculations,  which  are  made  at  a  single 
spatial  location  and  where  an  instrument  is  lowered  to  different  depths  to  calculate  an  attenuation 
profile. 


644  /  SPie  Vol  1302  Ocean  optics  X  (1990) 


From  a  visual  inspection  of  this  scatter  plot  for  «1,  it  appears  that  the  data  tend  toward  a 
horizontal  asymptote  as  z  increases  and  may  approach  a  vertical  asymptote  as  the  depth  nears 
zero,  although  having  no  data  for  depths  less  than  1.3  m  makes  this  assumption  subjective. 
Furthermore,  the  noise  and  spread  in  the  signals  of  LI,  L2,  and  L3  for  depths  less  than  3  m 
concomitantly  causes  noise  and  signal  spread  in  the  scatter  plots  of  /cl,  k2,  and  «3.  However, 
these  observations  and  the  fact  that  k  is  measured  in  inverse  meters  suggest  fitting  a  hyperbolic 
curve  of  the  form 


k  =  A/z  +  B 

through  the  data.  Moreover,  substituting  equation  (13)  into  equation  (6)  yields 


L  =  Lmin  +  exP(-2-(A/z  +  B)-(z-zr 


(13) 

(14) 


from  which  it  is  possible  to  solve  for  z  in  terms  of  L.  In  doing  so,  one  obtains  the  quadratic 
equation 


B-zz  +  (A  -  B-z  .  -  log(Y))-z  -  A*z  -  =  0, 
v  min  ov  "  min  ’ 


|  / 

where  we  have  let  Y  =  (AL/(L  —  Lmj  ))  ' 2  for  notations!  simplicity.  Using  the  quadratic 
formula,  we  obtain 

z  —  {l/2'B}*{log(Y)  +  B-zmjn-A.  [log(Y)-(l»g(Y)  +  2-D.zmin-2.A) 

+  +  A)*!7’}-  (16) 


There  are,  of  course,  two  solutions  to  the  quadratic  equation  expressed  in  (15V  It  turns  out  that 
in  taking  the  solution  corresponding  to  the  minus  sign,  we  obtain  negative  values  for  z  so  that 
this  solution  is  extraneous,  and  we  need  only  consider  the  solution  with  the  plus  sign. 

There  are  several  ways  to  estimate  the  unknowns  A  and  B.  One  way  would  be  to  select, 
presumably  judiciously,  two  points  (z,,^)  and  (z2,k2)  an^  determine  A  and  B  from 

simultaneously  solving  two  equations.  Alternatively,  one  can  estimate  B  by  looking  at  the  data 
for  depths  greater  than  20  m,  and  then  just  one  additional  point  (zpKj)  with  zx  small  is 

sufficient  to  calculate  A.  Another  way  would  be  to  perform  a  least-squares  fit  of  the  functional 
form  k  =  A/z  +  B  to  the  data. 

Yet  still  another  way  to  approach  the  determination  of  A  and  B,  which  produced  good 
results,  is  to  minimize  the  square  error  produced  from  equation  (15).  That  is,  given  a  data  point 
(zj.Lj)  from  one  of  the  bands,  substitution  of  (z^Lj)  into  the  quadratic  equation  (11)  will  yield, 

not  0,  but  some  presumably  small  error  Summing  the  squares  of  these  errors  we  have 


X[B'Zi2+(A-B-zmin-log(Yi)).zi-A.zmin]2=  £cf,  (17) 

where  again,  Yt  =  (AL/(Lj  -  Lmjn))1^2  for  notational  convenience.  By  taking  partial 
derivatives  with  repect  to  A  and  B  and  setting  them  equal  to  0,  we  obtain  the  equations 
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A'X(Zi-2min)J  +  B- 

A-£v(*i-W’  + 


-  W’  =  -  zmin>' 

B-  X’iM'i  “  zmin>2  =  Xl06(Y‘)'Z‘2  (Zi  “W’ 


(18) 

(19) 


which  are  easy  to  solve  for  A  and  B.  Applying  this  technique,  we  restrict  the  input  data  points 
(z^LJ  to  have,  for  example,  3  <  z{  <  10,  thus  avoiding  the  noisy  data  for  z{  <  3  and,  for 

Zj  >  10,  the  problem  of  small  variations  in  the  signal  Lj,  producing  large  variations  in  the  depth 

zv  The  upper  bound  may  be  changed  according  to  which  band  is  used;  we  used  10  for  LI  and  17 

for  L2  and  L3.  These  values  were  chosen  because  the  scatter  plots  of  LI,  L2,  and  L3  flatten  out  at 
approximately  these  depths;  thus,  the  greatest  variation  in  the  signal  versus  depth  occurs  for 
values  less  than  these  thresholds. 


In  Table  7,  results  of  employing  this  technique  are  given.  Figure  5  shows  the  curve 
generated  by  this  technique  plotted  over  the  data  for  bands  LI.  Figu^  6  shows  the  error  plotted 
versus  depth.  Note  that  for  water  depths  less  than  3  m,  the  error  is  large,  as  is  the  error  for 
depths  greater  than  10  m.  The  reason  for  this  is  the  noise  in  the  signal  in  the  shallowest  depths, 
and  for  the  deeper  water,  very  small  changes  in  the  signal  value  correspond  to  large  changes  in  the 
depth  value.  In  Figure  7  this  phenomenon  can  be  seen  as  well.  Similar  results  are  obtained  for 
bands  L2  and  L3,  with  the  error  becoming  larger  at  approximately  17  m. 


5.  CONCLUSIONS 

The  TUPS  data,  collected  in  situ,  display  the  expected  form  predicted  by  the 
Lambert-Beers  Law'.  Some  minor  variations  are  accounted  for  by  either  being  in  shallow  water 
(less  than  3  m)  where  the  signals  are  noisy,  or  in  deeper  water  (roughly  speaking,  greater  than 
10m)  which  appears  to  have  a  different  effective  attenuation  coefficient  than  that  for  water  depths 
in  the  3—  to  10-m  range  (the  deeper  water  is  clearer). 

In  an  attempt  to  allow  for  changes  in  the  effective  attenuation  coefficient  k,  we  have 
introduced  a  new  regression  model  based  on  k  having  a  hyperbolic  form.  This  technique  differs 
from  the  classical  single— band  bathymetric  regression  technique,  in  that  one  chooses  an  Lmax  at 

a  zmjn  value,  a  somewhat  subjective  choice.  Both  techniques  require  one  to  choose  a  value  for 

Lmin-  The  results  reported  here  show  that  both  methods  perform  equally  well  on  the  TUPS  data, 

with  the  new  technique  slightly  better  overall.  This  improvement  is  due  to  a  slightly  better  fit  in 
the  shallow  water  Hess  than  3  m)  and  because  the  bathymetry  curve  is  forced  to  pass  through  the 
point  (zmjn,  Lmax).  A  multiband  approach  can  likely  be  developed  using  this  new  regression 

technique.  This  development  may  improve  the  bathymetric  information  extraction  for  such 
noisier  data  as  LANDSAT  Thematic  Mapper  imagery. 
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Figure  1.2.  Normalized  upwelting  L2  *s  depth. 
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Figure  13  Normalized  upwellmg  L3  vs  depth 


Figure  2  1  Log  (LI)  vs  depth 
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Figure  3.  Empirical  radiance  curve. 
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Figure  4.  kl  vs  depth. 
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Figure  5.  LI  vs  depth. 
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Table  1.  TUPS  Sensor  Configuration,  June,  1988. 


CHANNEL 

SENSOR 

S/N 

1 

Temperature  sensor 

Sea  Bird,  Inc. 

SBE3-638 

2 

Conductivity  sensor 

Sea  Bird,  Inc. 

SBE4  234 

3 

465  nm  Upwelling  Light 
Biospherical  Instr. 

MCP-200H-7129 

4 

507  nm  Upwelling  Light 
Biospherical  Instr. 

MCP-200H-7130 

5 

532  nm  Upwelling  Light 
Biospherical  Instr. 

MCP-200H-7131 

6 

Transmissometer 

Sea  Tech,  Inc. 

165 

7 

Fluorometer,  Sea  Mar 
Tech,  Inc.  (signal) 

6000AR-235 

8 

Fluorometer,  Sea  Mar 
Tech,  Inc.  (scale) 

6000AR-235 

9 

488  nm  Downwelling 

Light,  Biospherical  Instr. 

QCP-200LM7101 

10 

Echo  Sounder 

Ulvertech,  LTD 

205 

Table  2.  TESS  Sensor  Configuration,  June,  1988. 


CHANNEL 

SENSOR 

1 

Voltage  Reference,  2  5V,  Analog  Devices 
AD580M 

2 

Pyroheliometer  1,  light  bulb,  Eppley  Instr. 
Model  50,  S/N  3039 

3 

Pyroheliometer  2,  hemisphere,  Eppley  Instr. 
Model  PSP,  S/N  8022D1 

4 

488  nm  Light  (linear),  Biospherical  Instr. 

Model  QCP-200HM-488,  S/N  7105 

5 

441  nm  Light  (log),  Biospherical  Instr.  Model 
QCP-200LM-441,  S/N  7102 

Table  3.  Correlation  Matrix  of  TUPS  and  TESS  Measurements. 


VARIABLE 

TESSPYROI 

TESSPYR02 

465  nm 

507  nm 

532  nm 

DEPTH 

TESSPYROI 

1 

■ 

TESSPYR02 

1 

465  nm 

r*— 

1HSE 

1 

I 

507  nm 

■■ 

0.989 

532  nm 

0.336 

■WSM 

0  975 

995 

1 

DEPTH 

0  195 

0  201 

-0.235 

B 

329 

-0.374 

1 

Table  4.  Single  Eand  Regression 
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Table  5.  Two  Band  Regression 


Table  6.  Three  Band  Regression 


LI,  L2,  L3 

DEPTH 

RM 

RMS 

0-5 

0.940 

0.3550 

3-5 

0.894 

0.2654 

0-10 

0.964 

0.6785 

3-10 

0.961 

0.6368 

0-15 

0.983 

0.6937 

3-15 

0.982 

0.6562 

0-17 

0.983 

0.7469 

3-17 

0.981 

0.7212 

0-22 

0.936 

2.297 

3-22 

0.928 

2.322 

Table  7.  Hyperbolic  k  Fits 
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